Using the acoustic peak to measure cosmological parameters 
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Recent measurements of the cosmic microwave background radiation by the Boomerang and 
Maxima experiments indicate that the universe is spatially flat. Here some simple back-of-the- 
envelope calculations are used to explain the result. The main result is a simple formula for the 
angular scale of the acoustic peak in terms of the standard cosmological parameters: t ~ 193[1 + 

§(i-no) + !(i-fc) + &fi A ]. 
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As we enter the era of precision cosmology, it gets in- 
creasingly difficult to understand how cosmological pa- 
rameters are extracted from observational data. The cos- 
mic microwave background radiation (CMBR) is a prime 
example. Anisotropics in the CMBR are influenced by a 
large number of cosmological parameters, including, but 
not limited to; the Hubble constant; the spatial curva- 
ture; the spatial topology; the vacuum energy density; 
the baryon density; the number of light neutrinos; and 
the amplitude and spectral index of the primordial den- 
sity perturbations. Using accurate maps of the CMBR, it 
should be possible to fix all of these parameters with great 
precision ffl |2| || . However, with so many parameters and 
so many physical effects to keep track of, it is hard to ex- 
plain how a particular parameter is extracted from the 
dataQ. The aim here is to provide a simple explana- 
tion of how the microwave background radiation can be 
used to measure the spatial curvature. The test was first 
proposed by by Doroshkevich, Zel'dovich and SunyaevQ, 
and has subsequently been developed by many authors. 
The most comprehensive treatment can be found in the 
work of Hu and White Q. 

The curvature measurement is based on simple geom- 
etry. If you know the physical size of an object and how 
far away it is, then by measuring its angular size you can 
infer the curvature of space. Suppose that the object has 
size A and is a distance B away. In flat space the angle 
subtended by the object is given by 



arccos 



A 2 \ 
l -2B 2 ) 



(1) 



However, if the space is negatively curved with radius of 
curvature R c , the angle will be given by 



= arccos 1 — 



cosh(A/R c ) ~ 1 
sinh 2 (B/i? c ) 



(2) 



Notice that (||) recovers that flat space result ([!]) in the 
limit R c — > oo. The expression for a positively curved 
space can be obtained by replacing R c by iR c in (||). 

To apply the angular size test we need to know the 
size of a distant object and how far away it is. When the 
test is applied to the microwave background radiation, 
the role of the distant "ruler stick" is played by the size 
of the sound horizon at last scatter, and the distance 
to the object is the radius of the last scattering surface. 
Both of these quantities depend on several cosmological 



parameters, but the spatial curvature turns out to be the 
dominant effect. 

We can measure the angular size subtended by the 
sound horizon by looking for a special feature in the 
CMBR angular power spectrum. The angular power 
spectrum is obtained by "fourier" analyzing the CMBR 
anisotropy pattern. Sound waves in the photon-baryon 
fluid with wavelengths roughly twice the size of the sound 
horizon at last scatter will have just reached a maximum 
density contrast when matter and radiation decouple. 
As we shall see, these waves have periods that are long 
compared to the time taken for matter and radiation to 
decouple. Thus, the compression-rarefaction pattern is 
snap frozen at decoupling. The enhanced density con- 
trast that occurs on the scale of the sound horizon leads 
to an enhanced temperature anisotropy, as the CMBR 
photons are redshifted by an amount proportional to the 
local density. By measuring the scale at which the peak 
in the angular power spectrum occurs, we are able to es- 
tablish the angular size of the sound horizon. Note: The 
acoustic peaks are density peaks, not "Doppler peaks" . 
The fluid is at a turn-around point when maximum den- 
sity contrast is reached, so the velocity of the baryons, 
and hence the Doppler shift, is at a minimum. 

Before proceeding to show how the angle subtended by 
the sound horizon is related to the spatial curvature, a 
little more non-euclidean geometry is in order. Consider 
a geodesic triangle drawn in hyperbolic space. The law 
of cosines reads 

cosh(C/i? c ) = cosh(A/R c ) cosh(B/R c ) 

- smh{A/R c ) sivh(B/R c ) cos 7 .(3) 

Here A, B, C are the side lengths and a, (3, 7 are the op- 
posite angles. Now suppose that C = B, B 3> A and 
R c S> A. Using the law of cosines we have 



(4) 



R c sinh(B/R c ) 



and 



= 



2 2i? c tanh( J B/i? c ) 
The sum of the angles in the triangle is given by 
E = a + 2/3 = 7T- a(cosh( J B/i? c ) - 1) + . . . . 



(5) 



(6) 
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Thus, the angle sum is less that 180° if space is negatively 
curved, greater than 180° if space is positively curved, 
and equal to 180° if space is flat. In our cosmological 
setting, A is the size of the sound horizon, B is the radius 
of the last scattering surface and a is the angular scale 
corresponding to the first acoustic peak. Space is flat if 
the angle sum in our cosmic triangle adds to 180°. 

Turning now from the spatial geometry to the space- 
time geometry, the unperturbed background geometry 
is described by the Friedman-Robertson- Walker line ele- 
ment 



ds 2 



-dt 2 + a(t) 2 (d X 2 + Rl sinh 2 (x/i? c )^ 2 ) . (7) 



Here a(t) is the scale factor in units where an = 1 today, 
and \R C \ is the spatial curvature radius today. The time- 
time component of Einstein's field equations reads (in 
units where G = c = 1) 



1 



8- 



i 2 R 2 



(8) 



where p is the energy density and a dot denotes dj dt. 
Solving for R c we find 



R,= 



1 



(9) 



where Hq — (d/a)o is the Hubble constant and f^o = 
Pa/pc is the total energy density today in units of the 
critical density p c — 3H 2 /8ir. Space is negatively curved 
if Slo < 1, positively curved if Hq > 1 and flat if £7q = 1. 
Using conservation of energy-momentum, the Friedman 
equation (j^) can be rewritten in the useful form 



i2 1( 



a 3(l+™) 



(10) 

Here fl r , £l m and Q w denote the contributions to the total 
energy density from radiation (photons and light neutri- 
nos) , non-relativistic matter (baryons and cold dark mat- 
ter), and an unclustered dark matter component with 
equation of state p — wp where w < —1/3. The un- 
clustered dark matter takes the form of a cosmological 
constant when w = — 1. 

Since angles are conformally invariant, we can use the 
conformally related static metric 



ds 



7P- 



-dr] 2 + d X 2 + R 2 C sinh 2 (x/i? c )^ 2 



(11) 



The static (or optical) metric has the nice property 
that null geodesies in spacetime correspond to ordinary 
geodesies in space. The conformal time rj is related to 
the cosmological time t by dr) = dt/a. Our task now is 
to calculate the size of the sound horizon at last scatter, 
and the radius of the last scattering surface. To lead- 
ing order, the sound speed in the photon-baryon fluid is 
equal to c s = 1 so the sound horizon is roughly 1 /v3 
times smaller than the conformal time interval between 



last scatter and the big bang (or between last scatter 
and reheating if we are considering inflationary models). 
Thus, the size of the sound horizon is given by 



Xsh 



1 

7i 



(Vsls - rjrh) 



1 



da 



(12) 



The size of the universe at last scatter, a s i si is inversely 
proportional to the redshift of last scatter, z s i s w 1100. 
The radius of the surface of last scatter is equal to the 
conformal time interval between last scatter and today: 



Xsls = V0 - Vsls 



da 
~Ha~ 2 



(13) 



Using the same approximations used to derive equation 
(0) , we can relate the angle subtended by the sound hori- 
zon, 9 sh , to Xsh and Xsls- 



Xsh 



R c sinh(xsis/Rc) 



(14) 



Let us begin with a simple case. Consider a matter 
dominated universe with Q r = £l w = 0. The integrals 
© and © yield 



Xsh — 



sis 



Xsh 



= R r arcsinh 



2yi - n 

fin 



(15) 



and the angular size of the sound horizon is given by 



sh — 




The angle sum in the cosmic triangle is given by 

« 1 - fio 



180° - 2° 



(16) 



(17) 



Note that while both Xsh and Xsls depend on Hq and 
fioj the angles only depend on fin- Thus, in a mat- 
ter dominated universe, the position of the first acous- 
tic peak is an excellent measure of the curvature. Since 
the CMBR experiments report their results in terms of 
angular power spectra, it is conventional to convert the 
angular scale into its fourier equivalent, the multipole 
number I ~ ir/9. For a matter dominated universe, the 
first acoustic peak in the angular power spectrum is lo- 
cated at £ peak w £ sh ~ 180/v^V 

For more realistic cosmological models, with both ra- 
diation and multi-component dark matter, the integrals 
( |l2| ) and (13) can not be evaluated in terms of simple 
functions. Approximate forms can be found as a power 
series expansions in the quantities fl w /Q mi il c /Q m and 
where J7 C = 1 — fio is a measure of the curva- 



Zsls I Zeq ; 

ture and 



O m /fi r denotes the redshift of matter- 



radiation equality. To leading order we have 



Xsh 



Zsls 



%sls 
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and 



Xsls 



(18) 



1 _ 1 Oc 

6 ri, r 



a, 



2(1 - 6w) ft r 



The expression for the radius of the last scattering sur- 
face, Xsls, is good to within 10% for |fi c /fi m | < 2 and 
|f2u,/f2 m | < 2 - see the appendix for details. The quantity 
z s is/ z eq that appears in ( |lS| ) is well approximated by 



1 



%sls 



(19) 



where h is the Hubble constant in units of 100 km s _1 
Mpc . In order to find a simple expansion for Xsh, we 
can choose either h or fl m h 2 as a free parameter and 
expand the quantity (y/l + z s i s /z eq - yjz s i s /z eq ) as 



1 

71 



(8n m h 2 - 1) 9(8n m h 2 - 1) 2 



2 / (i - n m ) 

3 V 10 



G4 



(1-h) 



The first version works best when Q m and h are close to 
the currently favored values of Q m = 0.3 and h = 0.65. 
The second version works best if Sl m ~ 1 and h ~ 1. 
Putting everything together in ([14]) we find 



1 



t _ in. 
2 a 



i 



, 2(1 - 6to) O r , 
(8^-1) , 



,(20) 



or, specializing to the case w = — 1 and fi TO ~ 1, we find 
V2 




1 



30 c 



(1 - h) 
35 5 



(21) 



The above expressions for 9 s h give a good qualitative 
picture of how the various cosmological parameters af- 
fect the location of the first acoustic peak. We see that 
the peak position is mainly determined by the curvature, 
and only weakly dependent on the value of the Hubble 
constant. The peak position is largely insensitive to the 
value of the cosmological constant. 

The angle subtended by the sound horizon is smaller 
in a negatively curved universe and larger in a positively 
curved universe. The angles in the triangle formed by 
the sound horizon and the Earth sum to 



E ~ 180°-1.15° 



1 



140, 



8n„,h 2 - 1 



(22) 

Neglecting photon self-gravity, sound waves in the 
photon-baryon fluid obey a simple harmonic oscillator 



equation, and the position of the first acoustic peak 
corresponds to the angular size of the sound horizon 
£ s h ~ n/Osh- However, when photon self-gravity is in- 
cluded, the oscillator equation gains an anharmonic term 
that shifts the position of the first few peaks. Taking this 
into account, and using standard isentropic initial con- 
ditions, and neglecting Silk damping, the temperature 
fluctuations vary as a function of scale as Q 

T w constant — cos^^) + — — sin(^^ s ^) . (23) 



IBsh 

Solving for the position of the first peak, we find 

7T 

**peak ~ 0.873 — , 



7sh 



so that for n r , 



1 



193 1 + 



35 



(1-h) 



(24) 



(25) 



5 35 5 

The recent Mat , Boomerang §] and Maxima^ re- 
sults locate the first acoustic peak at £ w 200, 1 — 197±6 
and I w 220, respectively. Using the Boomerang re- 
sults, and allowing h and Q\ to vary freely over the range 
0.5 < h < 0.8 and < fl\ < 0.8, our approximate for- 
mula ( p5[ ) yields a best fit value of f2g = 1-07 ± 0.1. 
This result is consistent with the universe being spa- 
tially flat, and agrees with the detailed analysis |l(| of 
the Boomerang data. Since the curvature is the dom- 
inant effect in fixing the location of the acoustic peak, 
we are able to get a good fix on Qq despite having only 
one equation for three unknowns (i/o, ^o> ^a)- I n con- 
clusion, simple analytic formulas can be found that give 
good qualitative, and decent quantitative, insight into 
how the CMBR observations are used to fix the spatial 
curvature. 
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appendix 

We made two major approximations in arriving at 
equation (|o|). The first was to treat the sound speed as 
a constant, when a more accurate approximation would 
be to set c s = 1/^/3(1 + £), where £ = 3af2b/4f2 7 is the 
baryon-photon momentum density ratio Keeping the 
next to leading term in £ s ; s , the size of the sound horizon 
is given by 



Xsh 



Csls 

"IT 



Zsis 



l 



. %sls 



Zsls 



ZsU 
/J p f. 



(26) 
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Assuming that there are three light neutrino species, the 
quantity £ s i s can be written as 



£,sls 




(27) 



For reasonable values of the cosmological parameters, we 
find £ s ; s < 1, so we can negelect the second term in equa- 
tion (fe 
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FIG. 1: The fractional error in the first order approximation 
to angular size distance to the surface of last scatter. The 
dashed lines mark contours of 8% error. The error is con- 
siderably less than 8% across most of parameter space. The 
missing corner corresponds to the portion of parameter space 
where Sl m = (Qo — ^a) < 0.2. 

The second major approximation was to expand the 



expression for r s i s = R c sinh(xsis / Rc) in terms of fl c /Q ri 
and fl w /il m : 



r s i s 



1 / fL 



a, 



8 



i a. 



2fl m 2(l~6w)n r 



1 , 



3 - 2w 



[i - i2w) \a 



4(1 - 2w)(i - 6w) \n m n 



(28) 



So long as fl m > Q c and tt m > fl w , the higher or- 
der terms can safely be neglected. Figure 1. shows 
the percentage error in the first order truncation of r s i s 
as compared to a full numerical evaluation. The frac- 
tional error is less than 8% across a wide portion of pa- 
rameter space, including the interesting region around 

(n c ,n A ) = (0,0.7). 

Our final task is to show that the period of the wave 
responsible for the first acoustic peak is large compared 
to the time taken for matter and radiation to decouple. If 
this were not the case, the anisotropy would not be frozen 
in and the acoustic peak would be washed out. The con- 
formal period of the wave is given by T ~ 2TTXsh/c s 
and the conformal time interval taken to decouple is 
A?y = rj{z s i s ) — r/(z s i s + Az), where Az w 300 is the 
redshift interval for decoupling. The ratio of T to A77 is 
given by 



Ar] \Az/y z eq y z eq \ z eq J J 

(29) 

For reasonable cosmological parameters, we find T ' j Ar\ > 
30, which tells us that the acoustic waves are effectively 
snap frozen when matter and radiation decouple. 
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